Critical properties of a dilute 0(n) model on the kagome lattice 
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S ; Abstract 

■ A critical dilute 0(n) model on the kagome lattice is investigated analytically and numerically. 

We employ a number of exact equivalences which, in a few steps, link the critical 0(n) spin model 
^ I on the kagome lattice to the exactly solvable critical g-state Potts model on the honeycomb lattice 

' with g = (n + 1)^. The intermediate steps involve the random-cluster model on the honeycomb 

o 
o 



lattice, and a fully packed loop model with loop weight n' = ^Jq and a dilute loop model with 
loop weight n, both on the kagome lattice. This mapping enables the determination of a branch 



J> ' of critical points of the dilute 0(n) model, as well as some of its critical properties. For n = 0, 

' this model reproduces the known universal properties of the Q point describing the collapse of a 

, polymer. For n 7^ it displays a line of multicritical points, with the same universal properties as 

' a branch of critical behavior that was found earlier in a dilute 0(n) model on the square lattice. 

00 ; 

■ These findings are supported by a finite-size-scaling analysis in combination with transfer-matrix 

•'-H . calculations. 
>< 

. ; PACS numbers: 05.50.+q, 64.60.Cn, 64.60.Fr, 75.10.Hk 
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I. INTRODUCTION 

The first exact results [l| for tlie 0(n) critical properties were obtained for a model on 
the honeycomb lattice, and revealed not only the critical point, but also some universal 
parameters of the critical state, as well as the low-temperature phase, as a function of n. 
The derivation of these results depends on a special choice of the 0(n) -symmetric interaction 
between the n-component spins of the 0{n) model, which enables a mapping on a loop gas 
2I. These results were supposed to apply to a whole universality class of 0(n)-symmetric 
models in two dimensions. 

Since then, also Oin) models on the square and triangular lattices were investigated 

n n 

pi, IJ] . Indeed, branches were found with the same universal properties as the honeycomb 
model, but in addition to these, several other branches of critical behavior were reported. 
Among these, we focus on 'branch 0' as reported in Refs. 4|. The points on this branch 
appear to describe a higher critical point. For = 0, it can be identified with the so- 
called 9 point describing the collapse of a polymer in two dimensions, which has been 
interpreted as a tricritical 0{n = model. It has indeed been found that the introduction 
of a sufficiently strong and suitably chosen attractive potential between the loop segments 
changes the ordinary 0{n = 0) transition into a first-order one sl, such that this change 
precisely coincides with the n = point of branch 0. Thus, the 6 point plays the role of 
a tricritical 0(?2 = 0) transition. Furthermore, it has been verified that tricriticality in the 
0(n) model can be introduced by adding a sufficient concentration of vacancies into the 

n 

system [7|]. More precisely, the introduction of vacancies leads to a branch of higher critical 
points, of which the points n = and n = 1 belong to universality classes (of the 6 point and 
the tricritical Blume-Capel model respectively) that have been described earlier as tricritical 
points. 

However, the critical points of branch on the square lattice appear to display universal 
properties that are different from those of the branch of higher critical points of the 0(?7,) 
model with vacancies [3], except at the intersection point of the two branches at n = 0. 
It thus appears that the continuation of the 6 point at n = to n 7^ can be done in 
different ways, leading to different universality classes. In order to gain further insight in 
this situation, the present work considers an 0{n) loop model on the kagome lattice with the 
purpose to find a 6'-like point, to continue this point to n 7^ and to explore the resulting 



universality. 



II. MAPPINGS 



The partition function of the spin representation of g-state Potts model on the honeycomb 
lattice 

^ Potts = 

^exp 

{S} \ <i,j> / 

depends on the temperature T by the coupling K = J/k-^T, where J is the nearest-neighbor 
spin-spin interaction. The spins Si can assume values 1, 2, ■ ■ ■ , g and their index i labels the 
sites of the honeycomb lattice. The first summation is over all possible spin configurations 
{S}, and the second one is over the nearest neighbor spin pairs. This Potts model can be 
subjected to a series of mappings which lead, via the random-cluster model and a fully- 
packed loop model, to a dilute 0(n) loop model which can also be interpreted as an 0(n) 
spin model. 



A. Honeycomb Potts model to fully-packed kagome loop model 

The introduction of bond variables, and a summation on the spin variables map the Potts 

n 

model onto the random-cluster (RC) model [8|, with partition function 

^RcKg) = $^«^'g^% (2) 

B 

where A*";, is the number of bonds, Nf. the number of clusters, and u = — 1 the weight of 
a bond. The sum is on all configurations B of bond variables: each bond variable is either 1 
(present) or (absent). In Eq. ([2]), q can be considered a continuous real number, playing 
the role of the weight of a cluster. Here, a cluster is either a single site or a group of sites 
connected together by bonds on the lattice. A typical configuration of the RC model on the 
honeycomb lattice is shown in Fig. [1] 

The next step is a mapping of the RC model on the honeycomb lattice onto a fully packed 
loop (FPL) model on the kagome lattice, which proceeds similarly as in the case of the square 
lattice j^. The sites of FPL model sit in the middle of the edges of the honeycomb lattice, 
and thus form a kagome lattice [lO| . Fully packed here means that all edges of the kagome 




FIG. 1: Mapping of the RC model onto a FPL model. The sites of the honeycomb lattice are 
shown as black circles. The dashed and the thin solid lines display the empty and the occupied 
edges (bonds) of the RC model on the honeycomb lattice respectively. The RC configuration is 
here represented by an FPL configuration on the surrounding lattice, i.e., the kagome lattice. Its 
loops (bold solid lines) follow the boundaries of the random clusters, both externally and internally. 
The Boltzmann weight of this finite-size configuration of the RC configuration is u^'^q^^ according 
to Eq. ([2]), and that of the corresponding FPL configuration is a^a^v?^ according to Eq. ([3]). 

lattice are covered by loop segments. The one-to-one correspondence between these two 
configurations is established by requiring that the loops do not intersect the occupied edges 
(bonds) of the honeycomb RC model, and always intersect the empty edges, as illustrated 
in Fig. [H 

To specify the Boltzmann weights of the FPL model, we assign a weight n to each loop, a 
weight ai to each vertex where the loop segments do not intersect an edge which is occupied 
by a bond of the RC model, and a weight a2 to each vertex where the loop segments intersect 
an edge which is empty in the RC model, as illustrated in Fig. [21 The partition function of 
the FPL model on the kagome lattice is thus defined as 

4^pi(ai, a2, n) = Y, «r^«rn"^' , (3) 
J" 

where mi is the number of type-1 vertices, m2 is the number of type-2 vertices and mi the 
number of loops. The sum is on all configurations T of loops covering all the edges of the 
kagome lattice. 





ai 



a2 



FIG. 2: Vertex weights of the FPL model. The bold solid lines represent loop segments. The 
weight of vertex where the loops do not intersect a bond (thin solid line) is oi. The weight of a 
vertex where two loops intersect an unoccupied edge (dashed line) is 02- 

The one-to-one correspondence between RC configurations and FPL configurations makes 
it possible to express the configuration parameters mi, m2 and mi of the FPL in those of 
the RC model, namely A^;, and Nc. Each vertex of type-1 corresponds with a bond of the 
RC model on the honeycomb lattice, thus 



mi = Nb. 



(4) 



The total number of the two kinds of vertices is equal to the number of edges on the 
honeycomb lattice, i.e.. 



mi + m2 



3N 



(5) 



where N is the total number of sites of the honeycomb lattice. Here we ignore surface effects 
of finite lattices. Furthermore, a loop on the kagome lattice is either one surrounding a 
random cluster on the honeycomb lattice, or one following the inside of a loop formed by 
the bonds of a random cluster. Thus 



mi = N, + Ni, (6) 
where Ni is the loop number of the RC model. Together with the Euler relation 

N, = N-Nb + Ni, (7) 
Eqs. (Hj) to ([6]) yield the numbers of vertices and loops on the kagome lattice as 

mi = A^;, 

m2 = 3N/2-Ni, (8) 
mi = 2Nc + Nh - N . 



Substitution in the partition function ([3]) leads to 



The weight of a given loop configuration is thus equal to the corresponding RC weight u'^'^q'^" 
if 

n = ^ 

ai = uq~^ (10) 



1 

02 = g3 , 



which completes the mapping of the RC onto the FPL model. 



B. Fully-packed loop model to dilute loop model 

Next we map the FPL model on the kagome lattice onto a dilute loop (DL) model on the 
same lattice, using a method due to Nienhuis (see e.g. Ref. |3|). The partition function of 
the FPL model on the kagome lattice is slightly rewritten as 

^^PL = («i + aa)"^ 5^<^<n(ri - 1) + ir (H) 
r 

with w\ = a\j (oi + a2) and W2 = (12/ (oi + 02). Eq. (ITT!) invites an interpretation in terms 
of colored loops, say red with loops of weight n — 1 and green loops of weight 1. Each of the 
2^1 terms in the expansion of [(n — 1) + 1]*"' thus specifies a way to color the loops: 

[(n - 1) + l]'"' = ^ {n-l)^n^\ 

{colorings} 

where Ir and Ig denote the number of red loops and green loops respectively, Ir + lg = mi. Let 
C denote a graph JF in which the colors of all loops are specified. The partition sum can thus 
be expressed in terms of a summation over all colored loop configurations C. The vertices 
of the kagome lattice are visited by two colored loops, and can thus be divided into 6 types, 
shown in Fig. [3] with their associated weights xi = yi = zi = wi and X2 = y2 = Z2 = W2- 
Thus, Eq. ( ITTj) assumes the form 

Zill = (a, + a^)- 5:xf-xf-yf-yf-.f-.f-(n - 1)'^'^ . (12) 
c 




Xl X2 Yi y2 Zl Z2 



FIG. 3: (color online). Weights of colored vertices. The vertical solid lines represent occupied 
edges (bonds) on the honeycomb lattice, while broken lines stand for empty edges. The bold solid 
lines represent the red loop segments, and the bold dashed lines the green ones. 

The sum on all colored loop configurations may now be replaced by two nested sums, 
the first of which is a sum on all dilute loop configurations of red loops, and the second 
sum J^gin is on all configurations of green loops Q that are consistent with TZ, i.e., the green 
loop configurations that cover all the kagome edges not covered by a red loop. Thus 

Z^A = («i + a,r-^J2-^'^-2-^^-^?-in - 1)'^ $:,f-yf-l'^ . (13) 
7^ g\n 

For each vertex visited by green loops only, there are precisely two possible local loop 
configurations. Since the loop weight of the green loops is 1, the summation over such pairs 
of configurations is trivial: 

E y^^'y^'' = E = (y^ + y^f' = i > (i^) 

g\n g\n 

where Ng is the number of green-only vertices. The FPL partition sum thus reduces to that 
of a dilute loop model, involving only red loops of weight n — 1: 

Zpp^l(«1'«2,^) = (oi + a2)^Z^^{xi,X2,Zi,Z2,n- 1) , (15) 

where the partition function of the dilute loop model is defined as 

Z^l{xi,X2,zi,Z2,n) = Ea;f"'a;f"''2;f"i2f''n^' , (16) 

c 

in which we forget the color variable, and denote the number of loops in a dilute configu- 
ration C a.s Ni. The dilute vertices are shown in Fig. |U together with their weights. The 
exponents of the vertex weights in Eq. ( fT6l) represent the numbers of the corresponding 




(a) (b) (c) 

FIG. 5: (color online). Partial summation on the green loops. The solid lines represent red loops, 
and the dashed lines green loops. For a fixed configuration of red loops, each vertex visited only by 
green loops has two possible weights: yi or 1/2 (see Fig. [3]). For the simple case shown here, there 
are two possible configurations (a) and (b), of which the relative weights are x\x2yiZiZ2{n — l)2l2 
and x\x2y2ZiZ2{n — 1)^1^ respectively. Addition of these weights yields the weight x\x2ZiZ2{n — 1)^ 
of the DL configuration shown in (c). 

vertices. Because of the similarity with the derivation of branch on the square lattice, we 
refer to the model flTBl) as branch of the kagome 0{n) loop model. 

The transformation between the FPL and the DL model is illustrated in Fig. [5l 

C. Dilute loop model to 0(n) spin model 

The Boltzmann weights in Eq. ( |T6|) contain, besides the loop weights, only local weights 
associated with the vertices of the kagome lattice. Just as in the case of the 0{n) model 
on the square lattice described in Ref. j^, there are precisely four incoming edges at each 
vertex. This implies that there is an equivalent 0(n) spin model: 

Z^l{Xi, X2, Zl, Z2, n) = Zspin(xi, X2, Zl, Z2) , (17) 
o 



of which the local weights have the same relation with the vertex weights as for the square 
lattice model of Ref. [Sj. Thus, the partition sum of the spin model is expressed by 

V 

)] • (18) 

The product is on all vertices v of the kagome lattice. The spins Sm sit on the midpoints 
of the edges of the kagome lattice. Their subscript "f 2" specifies the vertex v as well as the 
position i (with 1 < i < 4) with respect to the vertex. The label 1 runs clockwise around 
each vertex, such that the spins s^i and 8^2 sit on the same side of the honeycomb edge 
passing through vertex v. The spins have n Cartesian components and are normalized to 
length ^Jn. There are two different notations for each spin (because each spin is adjacent 
to two vertices), but a given subscript vi refers to only one spin. Here the number n is 
restricted to positive integers, of which only the case n = 1 is expected to be critical. 



^spin (^Xi, X2, ^1, Z2) 



D. Condition for criticality 



Since the critical point of the RC model on the honeycomb lattice is known ll|] as a 
function of g, namely 

(wy='-3g«)+g2 = 0, (19) 

the corresponding critical point of the n = ^/q FPL model on the kagome lattice is also 
known. According to Eq. ( fTOl) 

al = g5 , (20) 

from which the corresponding critical point of the DL model with loop weight n = — 1 
on the kagome lattice follows as 

™c _ _ ""he 



^2 = 4 = • (21) 



u 



he 



n 



III. DERIVATION OF SOME CRITICAL PROPERTIES 



The transformations described in Sec. [TTl leave (apart from a shift by a constant) the 
free energy unchanged, and lead to relations between the thermodynamic observables of the 
various models. Thus, the conformal anomaly and some of the critical exponents of the 
FPL and the DL models can be obtained from the existing results for the random-cluster 
■nodel. Thus. U.e in *e analogous case of the 0(„) n»de. on the s,ua.e lattice g the 
FPL model on the kagome lattice should be in the universality class of the low-temperature 
0{n) phase. However, the representation of magnetic correlations in our present cylindrical 
geometry leads to a complication. The kagome lattice structure, together with the FPL 
constraint, imposes the number of loop segments running along the cylinder to be even. 
Since the 0{n) spin-spin correlation function is represented by a single loop segment in the 
loop representation, which cannot be embedded in an FPL model on the kagome lattice, it 
is not clear how to represent magnetic correlations in this model. Thus we abstain from a 
further discussion of the scaling dimensions of the FPL model. 

1. conformal anomaly 

For the FPL model with loop weight n on the kagome lattice, the conformal anomaly c 



is equal to that of the n = ^/q Potts model 

mm-- 



1 — ^ -, 2 cos ^ = n, m>l. (22) 

m{m + 1) m + 1 



In the Coulomb gas language [M], it can be expressed as a function of the Coulomb gas 
coupling constant g, with g = m/ {m + 1): 

6(1 -o)^ 

c = l- 2cos{ng) = -n, 0<g<l. (23) 

The conformal anomaly c of the branch-0 critical 0{n) DL model on the kagome lattice 
with loop weight n is given by the same formula, but with n replaced by n -|- 1: 

c=l — ^ -, 2cos^^^ — = n + l, m>l. (24) 

m(m -|- 1) m+ 1 

The conformal anomaly is, via the number m, related to a set of scaling dimensions Xj 
as determined by the Kac formula [15|]: 

^ ^ [Pijm + 1) - qjm]'^ - 1 

2m(m + l) ' ^ ^ 



where Pi and are integers for unitary models. 

2. temperature exponent 

For the branch-0 critical DL model with loop weight n on the kagome lattice, the tem- 
perature exponent is expected to be the same as that for branch on the square lattice j^, 
namely Xf = Xj with pi = m,qi = m in Eq. fl25p . 

3. magnetic exponent 

The magnetic exponent of the branch-0 DL model with n = on the kagome lattice is 
not equal to the magnetic exponent of the low temperature 0(?2 + 1) loop model. The same 
situation was found earlier for the branch-0 0(0) model on the square lattice 3i. According 
to the reason given in [3|], the magnetic exponent is equal to the temperature one, i.e., the 
Pi = m,qi = m entry of Eq. fl2Sl) . The geometry of the underlying FPL model, where the 
number of dangling bonds is restricted to be even, plays here an essential role. Note that the 



a, 



even at the 9 point, is different 



magnetic exponent of the tri critical dilute 0(n) model 
from that of branch 0. 

These results for Xt and Xh are expressed in the Coulomb gas language as 

Xt = Xh = l- l/2g . (26) 



IV. NUMERICAL VERIFICATION 

A. Construction of the transfer matrix 



The transfer matrix is constructed for an L x M loop model wrapped on a cylinder, with 
its axis perpendicular to one of the lattice edge directions of the kagome lattice. The finite 
size L is defined such that the circumference of the cylinder is spanned by L/2 elementary 
hexagons (corner to corner). The cylinder is divided into M slices, of which L sites form a 
cyclical row, while each of the L/2 remaining sites forms an equilateral triangle with two of 
the sites of the cyclical row. The length of the cylinder is thus Ma/3. 

The partition function of this finite-size DL model is given by Eq. ([3]), but with Cm 
instead of C, in order to specify the length M of the cylinder: 

Z(M) ^ J]a;f-xf-^f-zf-n^' . (27) 

Cm 

1 1 



There are open boundaries at both ends of the cyhnder, so that there are L danghng edges 
connected to the vertices on row 1, as well as on row M . The way in which the end points 
of the dangling edges are pairwise connected by the loop configuration Cm is defined as the 
'connectivity', see Ref. Jsl for details. Here we ignore the dangling edges of row 1 (except 
for a topological property that will be considered later) and focus on the L dangling edges 
of row M . Since it is determined by the loop configuration, the connectivity (3 at row M is 
written as a function of Cm'- P = ({>{Cm)- The partition sum is divided into a number of 
restricted sums zj^\ each of which collects all terms in Z^'^^^ having connectivity (3 on row 
M, i.e.: 

2,„,^^2(M,_ 4«) = j:i,„„„,.f-.xf-^r-.^f-„''.. (28) 

13 Cm 

An increase of the system length M to M + 1 leads to a new configuration Cm+i which can 
be decomposed in Cm and the appended configuration Im+i on row M + 1. The graph Im+i 
fits the dangling edges of the loop graph Cm on the M-row lattice. The addition of the new 
row increases the number of the four kinds of vertices and of the number of loops by n^^^, 
71^2 J ^zii and rii respectively. The restricted partition sum of the system with M + 1 
rows is 



Cm+\ 

X]^ ^ ^ ^2 ' ^ ^ Sa,ifi{CM+i)-'^l ^-^2 ^^1 ^^2 ^ (29) 

Cm Im+i\Cm 

The last sum is on all sub-graphs Im+i that fit Cm- The connectivity ip{Cm+i) depends only 
on the connectivity /5 on row M, and on Im+i, so that we may write {p{Cm+i) = "ipiP, ^m+i)- 
Thus Eq. (1291) assumes the form 

Zif^^) = J2J2S,,,,CM)^f^^^^''4''4''n^' E S^,miM,.)<''^T'^T'4''n-'. (30) 

/3 Im+i\P 

The third sum depends only on a and (3, and thus defines the elements of the transfer matrix 
T as 

Tap= Yl ^»Mf3,iM+i)^T'x7'zi^'4''n'" , (31) 

Im+iW 

Substitution of T^/j and Eq. fl28|) in Eq. fl30|) then yields the recursion of the restricted 
partition sum as 

^(Am)^^y^^^(M)_ (32) 

/3 



In order to save memory and computer time, the transfer matrix of a system with finite 
size L is decomposed in ^ sparse matrices: 

T = Tl^^ ■ Tl^^_^ ■ Tl^^ ■ Tl ■ Tl^^ • . . . ■ T2 • Ti , (33) 

where Tj denotes an operation which adds a new vertex i on a new row, as illustrated in 
Fig. [61 Most of these sparse matrices are square, but Tlj^^ is not, because it increases the 
number of dangling bonds by two. The action of the other rectangular matrix, Tl^j^, reduces 
the number of dangling bonds again to L. 

During the actual calculations, we only store the positions and values of the non-zero 
elements of a sparse matrix, in a few one-dimensional arrays. Moreover, this need not be 
done for all the sparse matrices, because there are only four independent matrices. The 
other ones are related to these by the action of the translation operator j^, 

While the construction of the transfer matrix is formulated in terms of connectivities on 
the topmost rows M and M + 1, the connectivity on row 1 is not entirely negligible. In 
particular, the number of dangling loop segments on that row can be even or odd. As a 
consequence the number of dangling loop segments on the topmost row is then also even or 
odd respectively. This leads to a decomposition of the transfer matrix in an even and an 
odd sector. The odd sector corresponds with a single loop segment running in the length 
direction of the cylinder. 



B. Results of the numerical calculation 

For a model on an infinitely long cylinder with finite size L, the free energy per unit of 
area is determined by 

/W = ^lnAf, (34) 

where A^^'' is the largest eigenvalue of T in the = sector. From the finite-size data for 
f{L) we estimated the conformal anomaly c 



The magnetic correlation length ^h{L) is related to the magnetic gap in the eigenvalue 
spectrum of T as 

C'W = ^ln(Af/A?^), (35) 
where A^^"* is the largest eigenvalue of T in the rid = 1 sector. 



L-1 L 




(a) (b) 




(d) (c) 




(e) (f) 

FIG. 6: Constructing the transfer matrix. Appending a new row to the configuration is achieved 
in two parts. The first part consists of L/2 steps and is denoted • • - ^i (which are executed 
from right to left). Each step adds a new site to the lattice. Two of these steps are illustrated 
in (a) to (c). The number of dangling bonds does not change during these steps. The second 
part consists of L steps and is denoted ■ ■ •^L/2+i- The first step of these, Tl/2+i, adds a 

new vertex to the sub-row and increases the number of dangling bonds by 2 as shown in (d). The 
following steps " " ■^3L/2-i append vertices sequentially, and do not change the number of 

dangling bonds. After adding the last vertex by 13^/2 to the sub-row, the construction of a new 
row has been completed and the size of the system shrinks from L -|- 2 to L. 



The temperature correlation length C,t{L) is related to the temperature gap in the eigen- 
value spectrum of T as 

er^W = ^ln(A?VAf), (36) 

where A)^ is the second largest eigenvalue of T in the = sector. Using Cardy's con- 
formal mapping |l7|] of an infinite cylinder on the infinite plane, one can thus estimate the 
temperature dimension Xt and Xh- 

We calculated the finite-size data for the free energies of the FPL model at the critical 
points given by Eq. fl20l) for system sizes L = 2, 4, ■ ■ ■ , 28. These data include the case 
n = 0; this is possible because, for g — > one has u^^ = y/3q, so that the ratio between al 
and a2 in Eq. fl20l) remains well defined in this limit. 

The additional loop configurations allowed by the dilute model lead to a larger transfer 
matrix for a given system size, so that our results at the critical points given by Eqs. fl2T]) 
are restricted to sizes L = 2, 4, ■ • • , 18. The latter results also include the temperature 
and magnetic gaps. 

The finite-size data for the FPL and DL models displayed a good apparent convergence, 
and were fitted using methods explained earlier , [l6| , see also Ref . [l8| . 

In the kagome lattice FPL model, it is not possible to introduce one single open loop 
segment running in the length direction of the cylinder. The presence of a single chain 
would force unoccupied edges into the system, in violation of the FPL condition. Therefore, 
we have no results for Xh- Furthermore, in the case of the low-temperature 0{n) phase, 
the eigenvalue associated with Xt decreases rapidly when n becomes smaller than 2, and 
becomes dominated by other eigenvalues. Therefore, also results for Xf are absent for the 
FPL model, and our results are here restricted to the conformal anomaly c. The resulting 
estimates for the FPL model are listed in Tab. [H 

The results for the eigenvalue A^^'' of the the FPL model satisfy, within the numerical 
precision in the order of 10~^^, the relation between the FPL and DL models derived in 
Sec. IIIBI The larger dimensionality of the transfer matrix of the DL model in comparison 
with the FPL model generates new eigenvalues, and thus leads to new scaling dimensions 
that are absent in the FPL model. Final estimates for the conformal anomaly c and for 
the scaling dimensions Xf and Xh are listed in Tab. [TTl for the DL model. They agree 
well with the theoretical predictions, which are included in the table. Here we recall that, in 



TABLE I: Conformal anomaly c of the FPL model as determined by the transfer-matrix calculations 
described in the text. The sizes of the system L are from 2 to 28. Estimated error margins in 
the last decimal place are given in parentheses. The numerical results are indicated by 'num'. For 
comparison, we include theoretical values indicated by 'th', as given by Eq. ()23p . 
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1 
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analogy with the case of the branch-0 0{n) loop model on the square lattice [Sj, the magnetic 
scaling dimension should be exactly equal to the thermal one. This is in agreement with 
our numerical results. We found that the eigenvalues A^^'' and A^^"* were the same within 
the numerical error margin. Thus, we list only one column with results for the exponents in 
Tab. [Ill 



V. CONCLUSION 

We found a branch of critical points of the dilute loop model on the kagome lattice as a 
function of the loop weight n, which is related to the q = {n + l)^-state Potts model on the 
honeycomb lattice. The critical properties of these critical points are conjectured and verified 
by numerical transfer matrix calculations and a finite-size-scaling analysis. As expected, the 
model falls into the same universality class as branch of the 0(n) loop model Q| on the 

1 c 



TABLE II: Conformal anomaly c, magnetic scaling dimension and temperature scaling dimen- 
sion Xt of the DL model as determined by the transfer- matrix calculations described in the text. 
Estimated error margins in the last decimal place are given in parentheses. The numerical results 
are indicated by 'num'. For comparison, we include the theoretical values indicated by 'th', as 
given by Eqs. (IM]) and (f26D. 



n 




^num 






-1 


-2 


-2.0000 (5) 





0.0000000 (1) 


-0.75 


-1.3526699 


-1.3524 (3) 


0.073890718 


0.0738908 (2) 


-0.5 


-0.8197365 


-0.8194 (5) 


0.138570601 


0.138571 (1) 


-0.25 


-0.3749081 


-0.3747 (3) 


0.196602972 


0.196605 (5) 











1/4 


0.25000 (1) 


0.25 


0.31782377 


0.31778 (5) 


0.300602502 


0.30061 (5) 


^/2 - 1 


1/2 


0.500001 (1) 


1/3 


0.33334 (1) 


0.50 


0.58757194 


0.5876 (1) 


0.350604267 


0.35061 (1) 


\/3- 1 


4/5 


0.8002 (3) 


2/5 


0.3997 (5) 


0.75 


0.81497930 


0.8151 (3) 


0.404150985 


0.4037 (5) 


1 


1 


1.002 (3) 


1/2 


0.48 (3) 



square lattice. The analysis did, however, yield a difference. This is due to the geometry 
of the lattice. For the square lattice, it was found sl] that there exists a magnetic scaling 
dimension Xint,i as revealed by the free-energy difference between even and odd systems. 
Such an alternation is absent in the free-energy of the present model on the kagome lattice. 
While the number of dangling edges may be odd or even for the square lattice, it can only 
be even in the present case of the kagome lattice. 

The numerical accuracy of the results for the conformal anomaly and the exponents is 
much better than what can be typically achieved for an arbitrary critical point, whose loca- 
tion in the parameter space has to be determined in advance by so-called phenomenological 
I 

renormalization [19;]. This seems not only due to the limited precision of such a critical point. 
We suppose that the main reason is that irrelevant scaling fields tend to be suppressed in 
exactly solvable parameter subspaces. 



Acknowledgments 



We are much indebted to Bernard Nienhuis, for making his insight in the physics of 0(n) 
loop models available to us. This research is supported by the National Science Foundation 
of China under Grant #10675021, by the Beijing Normal University through a grant as well 
as support from its HSCC (High Performance Scientific Computing Center), and by the 
Lorentz Fund (The Netherlands). 



[1] B. Nienhuis, Phys. Rev. Lett. 49, 1062 (1982); J. Stat. Phys. 34, 731 (1984). 
[2] E. Domany, D. Mukamel, B. Nienhuis and A. Schwimmer, Nucl. Phys. B 190 [FS3], 279 
(1981). 

[3] H. W. J. B16tc and B. Nienhuis, J. Phys. A 22, 1415 (1989). 
[4] Y. M. M. Knops, H. W. J. Blote and B. Nienhuis, J. Phys. A 26, 495 (1993). 
[5] B. Duplantier and H. Saleur, Phys. Rev. Lett. 59, 539 (1987). 
[6] H. W. J. Blote, M. T. Batchelor and B. Nienhuis, Physica A 251, 95 (1988). 
[7] W.-A. Guo, B. Nienhuis and H. W. J. Blote, Int. J. Mod. Phys. C 10, 291 (1999); Phys. Rev. 
Lett. 96, 045704 (2006). 

[8] P. W. Kasteleyn and C. M. Fortuin, J. Phys. Soc. Jpn. 46 (Suppl.), 11 (1969); C. M. Fortuin 

and P. W. Kasteleyn, Physica (Amsterdam) 57, 536 (1972). 
[9] R. J. Baxter, S. B. Kelland, F. Y. Wu, J. Phys. A 9, 397 (1976). 

[10] I. Syozi, Progr. Theor. Phys. 6, 306 (1951). 
[11] F. Y. Wu, Rev. Mod. Phys. 54, 235 (1982). 

[12] H. W. J. Blote, J. L. Cardy and M. P. Nightingale, Phys. Rev. Lett. 56, 742 (1986). 
[13] I. Affleck, Phys. Rev. Lett. 56, 746 (1986). 

[14] B. Nienhuis, in Phase Transitions and Critical Phenomena, eds. C. Domb and J. L. Lebowitz 

(Academic Press, London, 1987), Vol. 11. 
[15] D. Friedan, Z. Qiu and S. Shenker, Phys. Rev. Lett. 52, 1575 (1984). 
[16] H. W. J. Blote and M. P. Nightingale, Physica A 112, 405 (1982). 
[17] J. L. Cardy, J. Phys. A 17, L385 (1984). 

[18] For reviews, see e.g. M. P. Nightingale in Finite-Size Scaling and Numerical Simulation of 



Statistical Systems, ed. V. Privman (World Scientific, Singapore 1990), and M. N. Barber in 
Phase Transitions and Critical Phenomena, eds. C. Domb and J. L. Lebowitz (Academic, New 
York 1983), Vol. 8. 
[19] M. P. Nightingale, Proc. Kon. Ned. Ak. Wet. B 82, 235 (1979). 



1 n 



